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Abstract 

We introduce a strategy to study the parameter space of the general, CP-conserving, two- 
Higgs-doublet Model (2HDM) with a softly broken ^-symmetry by means of a new “hybrid” 
basis. In this basis the input parameters are the measured values of the mass of the observed 
Standard Model (SM)-like Higgs boson and its coupling strength to vector boson pairs, the mass 
of the second CP-even Higgs boson, the ratio of neutral Higgs vacuum expectation values, and 
three additional dimensionless parameters. Using the hybrid basis, we present numerical scans 
of the 2HDM parameter space where we survey available parameter regions and analyze model 
constraints. From these results, we define a number of benchmark scenarios that capture different 
aspects of non-standard Higgs phenomenology that are of interest for future LHC Higgs searches. 
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1 Introduction 


The observation in 2012 m of a new boson with a mass close to 125 GeV [3j has been widely viewed 
as the discovery of the long sought after Higgs boson M- Despite the limited data set, the detailed 
analyses of the Run 1 LHC data by the ATLAS and CMS Collaborations have confirmed that the 
observed signal strengths (cross section times branching ratio) of the Higgs boson candidate relative 
to that expected of the Standard Model (SM) Higgs boson are consistent with SM predictions to 
within the accuracy of the current measurements w- Nevertheless, the current precision is at best 
20% in the bosonic channels and considerably less accurate in the fermionic channels. Thus, there is 
still considerable room for new physics beyond the SM to contribute to the properties of the already 
discovered scalar state. 

The Standard Model posits that the dynamics of one complex hypercharge-one, weak doublet of 
scalar fields is solely responsible for electroweak symmetry breaking. Three of the four degrees of 
freedom of this scalar doublet provide the longitudinal modes of the massive W± and Z gauge bosons. 
The remaining scalar degree of freedom is the SM Higgs boson. But, why should the scalar sector of the 
SM be of minimal form? The spin-1/2 quark and lepton degrees of freedom of the SM appear in three 
generations. The origin of the non-trivial flavor structure of the SM is presently unknown. By analogy, 
one might expect a replication in the scalar sector as well—a non-minimal Higgs sector consisting of 
multiple doublets. Adding additional doublets yields new phenomena in the scalar sector—charged 
Higgs scalars, neutral Higgs scalars of opposite CP quantum numbers (if CP is conserved by the scalar 
potential) or neutral Higgs of indefinite CP (if CP is violated by the scalar potential). The current 
Higgs data do not rule out the existence of an extended Higgs sector which, as we will discuss in some 
detail, could be realized at mass scales close to (or in the extreme case even below) 125 GeV. Thus, 
apart from any other theoretical motivation, it behooves us to devote a dedicated program at the LHC 
to search for evidence of a non-minimal Higgs structure. 

Of course, theoretical arguments have been also advanced in support of a non-minimal Higgs 
sector. Perhaps the most persuasive is based on the fact that the SM is unnatural mm -namely, it 
is difficult to understand how the scale of electroweak symmetry breaking [y = 246 GeV) arises in a 
more fundamental framework that includes gravity and its associated Planck scale, Mpl — 10 19 GeV. 
In the context of the SM coupled to gravity, electroweak symmetry is achieved only by fine-tuning 
the squared-mass parameter of the scalar potential to an accuracy of 34 decimal places. New physics 
beyond the SM that attempts to address this issue introduces new phenomena that enters at or near 
the TeV scale. Many such approaches invoke non-minimal Higgs sectors. The most well studied 
example of this kind is the minimal supersymmetric extension of the SM (MSSM), which employs a 
two-Higgs doublet scalar sector with quartic terms in the scalar potential that respect supersymmetry 

(susy) mm - 

Since we do not presently know the precise nature of the new physics beyond the SM, it is prudent to 
be open minded about the possible structure of the non-minimal Higgs sector. However, the observed 
value of the electroweak parameter p = myy/{m 1 2 z cosyy) ~ 1 strongly favors non-minimal Higgs sectors 
comprised only of singlet and doublet scalar fields }15)ll6j PI Among theories in this category, we find 
the two-Higgs doublet modelEI (2HDM) particularly attractive as the minimal extension that yields 
new scalar phenomena, including new charged states (H±) and neutral states with different (or mixed) 

1 Including triplet scalar fields for example typically violates p ~ i pmh] except in special cases that must be 
considered as fine-tuned m- 

2 For a comprehensive review of the 2HDM, see ng. A review that treats the 2HDM in the formalism employed in 
this paper can be found in m- 
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CV properties. Direct searches for the additional 2HDM Higgs states based on run-I LHC data have 
appeared both from ATLAS [21] and CMS [221123] . Interpretations of the recently discovered Higgs 
boson at 125 GeV in the context of the 2HDM and implications for future LHC searches have been 
presented in I24lf34] . 

In this paper, we develop a set of benchmark scenarios for LHC Higgs searches that capture 
different aspects of 2HDM phenomenology. Our scenarios are devised taking into account that the 
observed Higgs boson already possess properties that are close to those expected in the SM. The 
experimental data already provide important constraints on the 2HDM framework. For example, the 
absence of observable flavor changing neutral currents (FCNCs) mediated by tree-level exchange of 
neutral Higgs bosons requires that the structure of the Higgs-fermion Yukawa couplings must be of a 
special form [35, 36]. The simplest way to eliminate tree-level Higgs-mediated FCNCs is to impose a 
particular discrete Z 2 symmetry [cf. Table [Din Section \2. 4] that is exactly respected by the dimension- 
four quartic terms of the Higgs Lagrangian |371l38 j FI We make one further assumption in our analysis 
by imposing CV symmetry on the scalar potential. This assumption is not required based on data 
(even though a pure CV -odd state has been ruled out as explanation for the 125 GeV state [39]). The 
restriction to CP-symmetry provides an additional simplification to the analysis. In particular, the 
neutral scalar spectrum consists of two CP-even states h and H with mh < mn, and the CP-odd 
state A. Indeed, the methods developed in this paper are rather easily generalized to a 2HDM with a 
CP-violating scalar potential or CP-violating vacuum. This possibility will be addressed elsewhere. 

In Section 2, we provide the theoretical background relevant to the 2HDM. We introduce the Higgs 
basis mu. which is especially useful in our analysis, as it provides a very clear way of parametrizing 
the terms that yield deviations of the properties of the observed Higgs boson from those of the SM. 
Ultimately, we construct a “hybrid basis” of parameters, which we will employ in 2HDM parameter 
scans. These parameters include the two CP-even neutral Higgs boson masses (m/j and mn, where one 
of these masses is identified with the mass of the observed Higgs boson), the parameter cg_ Q which 
parametrizes the deviation of the lightest CP-even neutral Higgs boson h from that of the SM, the 
ratio of neutral Higgs vacuum expectation values, tan /3, and three dimensionless quartic couplings of 
the Higgs basis. Two of these three dimensionless couplings are related in a simple way to the masses 
of the charged Higgs and neutral CP-odd Higgs boson of the 2HDM. 

In Section 3, we present our numerical scans of the 2HDM parameter space using the hybrid basis 
of parameters. Based on these scans, which are performed using the code 2HDMC [421143] . we develop 
seven different benchmark scenarios: (1) a SM-like h, with cp_ a small but non-zero as allowed by 
current data to yield interesting phenomenology of the heavier H\ (2) a “flipped” scenario with a 
SM-like H (which implies that |cg_ Q | is near 1); (3) overlapping CP-even and CP-odd scalars with 
masses around 125 GeV; (4) a SM-like h and a heavy non-SM-Higgs mass spectrum with short cascade 
decays H —» ZA or H —>• (5) a SM-like h and a heavy non-SM Higgs mass spectrum with 

long cascade decays P± -> W ± A -»• W ± ZH or A -»• -> W^W^H-, (6) SM-like hVV and hit 

couplings (V = and Z), but with opposite sign hbb and hr + T~ couplings as compared with the 
SM; and (7) an MSSM-like scenario, in which the scalar potential parameters (with one exception) are 
fixed by the corresponding tree-level MSSM relations. Finally, in Section 4, we present our conclusions. 
The proposed benchmark scenarios are summarized in a set of tables presented in Appendix A. 

3 Dimension-two quadratic terms in the Higgs Lagrangian that softly break the Z 2 discrete symmetry are allowed since 
these terms will only generate Higgs-mediated FCNCs at the loop level, which typically are not large enough to be in 
conflict with experimental data. 
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2 Theoretical background 


2.1 The general two-Higgs doublet model (2HDM) 

The two-Higgs doublet model (2HDM) consists of two complex SU(2) doublet, hypercharge-one fields, 
4>i and !> 2 , and an SU(2)xU(l)-invariant scalar potential, 

V = +m| 2 <t 4$ 2 - [m \ 2 §[$2 + h.c.] + iAi^i ) 2 + ±A 2 (^$ 2 ) 2 + A 3 (${$i)($£$ 2 ) 

+A 4 ($I$ 2 )($5$i) + {|A 5 (^I $ 2 ) 2 + [A 6 ($ 1 $i) + A 7 ($£$ 2 )]4>1 $ 2 + h.c.} , ( 1 ) 

where the parameters rnj 2 , X$, Xq and A 7 are potentially complex. All other scalar potential parameters 
are manifestly real. We assume that the parameters are chosen such that the minimum of the scalar 
potential spontaneously breaks the SU( 2 )xU(l) electroweak gauge symmetry to U(1)em jSMS]. That 
is, at the minimum of the scalar potential, the neutral components of the complex doublet scalar fields 
acquire vacuum expectation values (vevs), 

= (i = 1,2), (2) 

where by convention v\ and u 2 are real and non-negative. The combination v 2 = v 2 + v 2 — (246 GeV ) 2 
is fixed by its relation to the Fermi constant and the W boson mass, v 2 = l/(\/2 Gf) = 

In the most general 2HDM, the fields 4*1 and d >2 are indistinguishable. Thus, it is always possible 
to define two orthonormal linear combinations of the two scalar doublet fields without modifying any 
prediction of the model. Performing such a redefinition of fields (henceforth called a change of basis 
of the scalar doublet fields ) leads to a new scalar potential with the same form as Eq. (JT|) but with 
modified coefficients. In this paper we shall focus on the case where the scalar potential and the 
vacuum state are CP-conserving, leaving the more general case for future work. 

The scalar potential is explicitly CP-conserving if and only if there exists a basis choice for the 
scalar fields in which m 2 2 , A 5 , A 6 and A 7 are simultaneously real. The general conditions that guarantee 
the existence of such a basis (called a real basis ) were developed in [56]. Henceforth, we shall assume 
that all scalar potential parameters given in Eq. (JT]) are real. However, it is still possible that the 
vacuum spontaneously breaks CP. In particular, spontaneous CP-violation takes place if and only if 
the scalar potential is explicitly CP-conserving, but no real basis exists in which the scalar vacuum 
expectation values are simultaneously real. Sufficient conditions for a CP-conserving vacuum are easily 
obtained. The minimization of the scalar potential fixes v\ = v cos /3, u 2 = v sin/3 (where 0 < /3 < ^ 7 r) 
and the relative phase of the two vevs, £ = £2 — £1 through the equations, 


mh = m\ 2 tan /3 cos £ — ^v 2 [Aicjg + (A 3 + A 4 + A 5 cos 2 £)s| + 3A §spcp cos £ + A 7 tan /3 cos£], (3) 

m 2 2 = iti 2 2 cot /3 cos£ — \v 2 [\+ (A 3 + A 4 + A 5 cos2£)pg + cot /3 cos£ + 3X7 spcp cos£], (4) 

ml 2 sin £ = \v 2 [2A 5 spcp cos £ + A 6 c^ + A 7 s^] sin £ , (5) 

where sp = sin/3 and cp = cos/3 and all scalar potential parameters are real by assumption. If 
sin£ 0 , then Eq. ([5]) can be used to obtain 


cos£ = 


m i 2 ~ ~ ^7 v 2 

X5VW2 


( 6 ) 
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Moreover, this value of £ corresponds to the minimum [maximum] of the scalar potential if A 5 > 0 
[A 5 < 0]. Thus, it follows (e.g., see Appendix B of ref. [37]) that no scalar potential minimum occurs 
for sin £ 7 ^ 0 if 

\m 2 l2 - ~ 2 > A 5 U 1 U 2 , (7) 

Thus if Eq. (0) is satisfied, then the scalar potential is minimized for sin£ = 0 (i.e., £ = mr for 
integer n). One can then perform a hypercharge gauge transformation, <f>j —>■ 4>,;, followed by the 

field redefinition, 4>2 —> (—l) n< I> 2 , which yields real non-negative vevs = Vi/y /2 (for i = 1 , 2 ). 

If Eq. (0) is not satisfied, then the scalar potential is minimized for sin£ 7 ^ 0. Nevertheless, one 
still must check that there exists no basis transformation of the scalar fields such that the resulting 
vevs are real. This can be accomplished without explicitly considering all possible basis choices for the 
scalar fields by making use of the so-called Higgs basis. In Section 12.21 we define the Higgs basis and 
explain how it can be used to verify that a 2HDM with an explicitly CP-conserving scalar potential 
also has a CP-conserving vacuum. 

2.2 The Higgs basis 

It is convenient to define new Higgs doublet fields, 

(Ht\ _ vie~ l ^ l <bi + u 2 e - ^ 2 <3?2 TT (H$\ _ -v 2 e l &$i + uie^ 1 d > 2 , oN 

= (//?) =- v -’ 2= {h?J= - v -' (8> 

It follows that (Hf) = v/y/2 and (H 2 ) = 0. This is the Higgs basis [40IHT ], which is uniquely defined 
up to an overall rephasing, H 2 —>• e lx H 2 . 

In the Higgs basis, the scalar potential takes the same form as Eq. (jTJ) but with new coefficients, 


V = Y\H\h x + Y 2 h\h 2 + [Y 3 H{H 2 + h.c.] + \Z x [h\h x ) 2 

+\z 2 (hIh2) 2 + z 3 (h\h 1 )(hIh 2 ) + z a (h\h 2 )(hIh x ) 

+ {\Z b (H\H 2 f + [Zq(h\h{) + Z 7 {h\H 2 )\h\h 2 + h.c.} , (9) 

where Y\. Y 2 and Z A ,. .., Z 4 are real and uniquely defined, whereas I 3 , Z 5 , Zg and Z 7 are complex 
and transform under the rephasing of H 2 , 

[Y 3 ,Z 6 ,Z 7 \^e- ix [Y 3 ,Z 6 ,Z 7 \ and Z 5 ->■ e~ 2ix Z 5 . (10) 

For an explicitly CP-conserving scalar potential where all the coefficients given in Eq. © are real, the 
real coefficients of the scalar potential in the Higgs basis are given by 

Y\ = + m\ 2 s 2 p — m\ 2 s 2 g cos £ , ( 11 ) 

Y 2 = m 2 u s 2 p + rn 2 2 Cp + m\ 2 s 2 g cos £ , ( 12 ) 

Z\ = \\Cp + A 2 s'jj + ^(A 3 + A 4 + A 5 cos 2 £)s 2/3 T 2 s 2 / 3 (^ 6 c /3 T A 7 s 2 p) cos £ , (13) 

Z 2 = Ais^ + A 2 C^ + ^(A 3 + A 4 + A 5 cos 2£)s^ ~~ 2 s 2 / 3 (Agsjg + \ 7 (?g) cos £ , (14) 

Z 3 = |(Ai + A 2 — 2 A 3 — 2 A 4 — 2 A 5 cos 2 £)s 2 j a + A 3 — (Ag — \ 7 )s 2 gc 2 g cos£ , (15) 

Z 4 = |(Ai + A 2 — 2 A 3 — 2 A 4 — 2 A 5 cos 2 £)s 2/3 + A 4 — (Ag — \ 7 )s 2 / 3 C 2 g cos£ , (16) 


5 




where s 2 p = sin 2/3, c 2 p = cos 2/3, etc. The potentially complex coefficients of the scalar potential in 
the Higgs basis are given by 

Y 3 = -e~^ [\irnh - m% 2 )s 2/ 3 + m 2 2 c 2(3 cos£ + im\ 2 sin£] , (17) 

Z 5 = e -2 ^ j \ [Ai + A 2 - 2(A 3 + A 4 + A 5 cos 2£)] s 2/3 + A 5 cos 2£ - (A 6 - A 7 )s 2/3 c 2/3 cos £ 

+i [A 5 c 2/3 sin 2 £ - (A 6 - A 7 )s 2( g sin £] \ , (18) 


Zq = e ^ j — ^[Aic^ — \ 2 s 2 p - (A3 + A4 + A5Cos2£)c 2j g]s 2( 3 + (\qC/3C 3 /3 + A 7 s ( gs3 j a)cos^ 


+*[5A 5 s 2/3 sin 2£ + (A 6 c| + A 7 s|)sin^] 


(19) 


Zj = e ^ j — \[\is 2 p — A 2 c 2 p + (A3 + A4 + A5COS2£)c 2j g]s 2j g + (A6S^S3/3 + A 7 c i ac3 ( a)cos^ 

+i[-lA 5 s 2 / 3 sin2£ + (A 6 s^ + A 7 c|)sin£] | , (20) 

One can check that by using Eqs. Q Q, one recovers the expected scalar potential minimization 
conditions in the Higgs basis, 

Y 1 = -±Z 1V 2 , Y 3 = -\Zqv 2 . (21) 

Note that if sin£ = 0 then Y 3 , Zq, Zq and Z 7 are all real and the scalar potential and the vacuum 
are CP-conserving. More generally, if all the coefficients of the Higgs basis are real for some choice 
of x [cf. Eq. (HOI) ] , then it follows that any basis related to this real Higgs basis by a real orthogonal 
transformation of the two doublet fields is also a real basis with real vevs. Thus, it follows that the 
scalar potential and vacuum are CP-conserving if and only if 11 TTi) 17b1T)l 

Im (ZtZl) = Im(Z 5 *Z 2 ) = Im(Z*Z 7 ) = 0 . (22) 

No separate condition involving Y 3 is needed in light of Eq. (ED- Thus, if sin^ 7 ^ 0, then it may 
still be possible that the vacuum is CP-conserving if Eq. (1221) is satisfied. One can obtain the general 
conditions for a CP-conserving vacuum by inserting Eqs. (1181) (1201) into Eq. (1221) . For example, 

Im(ZgZ 7 ) = j sin2^s 2 ^ [As(Ai—A 2 )—Aq+A 7 ] —^ sin^s 2 ( gc 2j a [AiA 7 +A 2 Ag —(A 3 +A 4 +A 5 )(A 6 +A 7 )] . (23) 

One noteworthy observation concerns the special case of Ai = A 2 and A 7 = — Xq (dubbed an exceptional 
point of the 2HDM parameter space in Ref. |3T]). I n fact, if both these relations hold simultaneously 
in one basis then they hold simultaneously in all bases. Applying this observation to the Higgs basis, 
it follows that Z\ = Z 2 and Z 7 = —Zq, which can be easily verified using eqs. m , (ED, ED and 
m■ As a result Im(ZgZ 7 ) = —Im(|Zg| 2 ) = 0, which is confirmed by Eq. (1231) . 

The general expressions for Im(ZgZ|) and hn(Z^Zj) are much more complicated and not partic¬ 
ularly illuminating. Nevertheless, if Eq. (|22l) is satisfied for some value of £ 7 ^ 0, then it follows that a 
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real Higgs basis exists. Performing an 0(2) basis transformation then yields a scalar potential of the 
form given by Eq. ([!]) with real vevs. That is, without loss of generality, we can assume that sin£ = 0 
in Eqs. (ITT]) (l20|) . which yields a real Higgs basis that is unique up to a possible sign ambiguity that 
corresponds to redefining the Higgs basis field H 2 —>■ —H 2 . 

Thus, we shall parametrize the general CP-invariant 2HDM scalar potential by its real Higgs basis 
form given in Eq. © , where all the Y) and Z* are real. The scalar potential minimum conditions given 
in Eq. (l2ljl fix Y 4 and Y 3 in terms of Z\ and Zg, respectively. Since (H±) = v/V2 and (H®) = 0, it 
follows that Y 2 , Z\, ..., Z 5 and the product Z^Z-j are uniquely defined, whereas Z§ and Z 7 separately 
change sign under the only possible transformation among real Higgs bases, H 2 —>• — H 2 . That is, Z$ 
and Z 7 are pseudo-invariant quantities with respect to the real Higgs basis transformation H 2 —> —H 2 . 
However, the relative sign of Z§ and Z 7 is meaningful. 

Physical observables must be invariant with respect to any possible Higgs basis transformation. 
In the most general 2HDM, tan/3 is also unphysical [411150] since there is no physical significance to 
an arbitrary real basis of scalar fields (apart from the Higgs basis and the basis defined in terms of 
neutral Higgs mass eigenstates). However, it is often the case that the form of the Higgs-fermion 
Yukawa couplings pick out a special scalar held basis, in which case tan/3 is promoted to a physical 
parameter. 

It is convenient to employ the real Higgs basis to evaluate the spectrum of Higgs masses in the 
CP-conserving model. The physical charged Higgs boson is the charged component of the Higgs-basis 
doublet H 2 , and its mass is given by 

m 2 H± = Y 2 + \Z 3 v 2 . (24) 

The three physical neutral Higgs boson mass-eigenstates are determined by diagonalizing a 3 x 3 real 
symmetric squared-mass matrix that is defined in the Higgs basis [401150] 

/ Z\v 2 Zqv 2 0 \ 

M 2 =[Z 6 v 2 Y 2 + ±(Z 3 + Z 4 + Z 5 )v 2 0 . (25) 

V 0 0 Y 2 + ±(Z 3 + Z 4 -Z 5 )v 2 J 

We immediately identify the CP-odd Higgs boson A = \/2lm H$ with squared mass, 

m\ = Y 2 + \{Z 3 + Z 4 — Z 3 )v 2 . (26) 


Note that the real Higgs mass-eigenstate held A is defined up to an overall sign change, which corre¬ 
sponds to the freedom to redefine H 2 —> — 3 ^ 2 - 

The upper 2x2 matrix block given in Eq. (1251) is the CP-even Higgs squared-mass matrix, 


M 2 h = 


Z\v 2 

Z 6 v 2 


Z 6 v 2 

m\ + Z 5 v 2 


(27) 


where we have used Eq. (1261) to eliminate Y 2 . 
eigenstates, h and H (with rtih < mjj) by 


To diagonalize M 2 H: we define the CP-even mass 

(28) 


(H\ = ( Cfl—a -s t 3 -cA A/2 Re Hi - v\ 

V h) c 8-a) V V^Re H$ ) 


where cp- a = cos(/3 — a) and sp- a = sin(/3 — a) are dehned in terms of the mixing angle a that 
diagonalizes the CP-even Higgs squared-mass matrix when expressed in the original basis of scalar 
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fields, {<I>i, $ 2 }- The real Higgs mass-eigenstate fields H and h are defined up to an overall sign 
change. This implies that f3 — a is defined modulo it. In particular, 


^ C/3—a > 

S/3—a ^ $(3—a 5 

h 2 -a +h 2 

=> H^-H, 

h —y — h , 

(29) 

—)• +Cg_ Q , 

S/3—a ^ $(3—a 5 

h 2 -A -h 2 

=> H^+H, 

h —y — h , 

(30) 

—Cp-a j 

s (3-a + 5 /3-a 5 

h 2 -A —H 2 

=> H -A —H , 

h —y -\-h . 

(31) 


It follows that the product sp- a cp- a is a pseudo-invariant quantity with respect to the real Higgs 
basis transformation H 2 —>• — H 2 . 

The squared masses of h and H are then given by, 

m H,h. = \ + ( z i + z 5 )v 2 ± \J \m 2 A + (Z 5 - Zi)v 2 ] ~ + 4 Z%v 4 | . (32) 


The following identity therefore holds, 

\Zq\v 2 = (m 2 H - Zxv^Zxv 2 - ml) . (33) 

Hence, diagonalizing M 2 H yields the following expressions: 

Ziu 2 = m 2 h s 2 p_ a + m 2 H c 2 p_ Q , (34) 

Zqv 2 = (m\ - m 2 H )sp- a c/ 3 - a , (35) 

2 1 ry 2 _ 2 2 1 2 2 (oc\ 

m A + Z 5 V — r nH s /3-a + m h c /3-a > (36) 

where sp_ a = sin(/3 — a) and cg_ Q = cos({3 — a). Note that Eq. ([35]) imply that 

Zq sp-a.cp-a < 0 . (37) 


Indeed, Eq. (1571) is invariant with respect to the real Higgs basis transformation H 2 -A- —H 2 , so the 
sign of the quantity Z§sp_ a cp_ a is physically meaningful. 

Using the fact that /3 — a is defined modulo it, we shall establish a convention where 


0 < /3 — a < tt . 


(38) 


In this convention, a is non-negative, the sign of the field h is fixed and cs~ a is pseudo-invariant 
with respect to the real Higgs basis transformation H 2 —>■ —H 2 . One can then derive expressions for 
cg_ Q and sp_ a from Eqs. (IM1) and ([35]) . where the signs of the corresponding quantities are fixed by 
Eqs. (1371) and ([38]) : 


C/3 —a 


sgn (Zq) 


/ 


Z\v 2 — m\ 
m H~ m h 


_ ~Z 6 v 2 _ 

\ m H ~ m h)( m H ~ Z lV 2 ) 


(39) 


S(3—a — 


/ 


m? H — Z\v 2 
m H~ m h 



(40) 


Note that we have used Eq. (1331) to obtain the second form for cg_ Q and sp_ Q in Eqs. (1391) and (1401) . 
respectively. 












2.3 The CP-conserving 2HDM scalar potential with a softly-broken discrete Z 2 
symmetry 


When we introduce the Higgs-fermion Yukawa couplings in Section 12.41 there will be some motivation 
to restrict the parameter freedom of the most general scalar potential given in Eq. (0) by requiring 
the invariance of the scalar potential under the discrete Z 2 symmetry $1 —>• +3>i and <f > 2 —7 — $ 2 . 
Imposing this discrete symmetry implies that m ? 2 = Ag = A 7 = 0 in Eq. 0). In this case, A 5 is the 
only potentially complex scalar potential parameter, which can be rendered real by an appropriate 
rephasing of $ 4 . It then follows from Eqs. Q and 0 that the scalar potential is minimized for 
sin£ = 0 if A 5 < 0 and cos£ = 0 if A 5 > 0. In the latter case, ($ 2 )/($i) = ±*tan/3. However, a 
redefinition of <3?i —7 yields real vevs while A 5 —7 — A 5 . Thus, the Z 2 -invariant scalar potential 

and the vacuum are CP-invariant. 

One can relax the discrete symmetry by allowing for m ? 2 7 ^ 0 in Eq. (JTJ) which softly breaks the Z 2 
symmetry. The quartic terms in Eq. (JTJ) still respect the Z 2 symmetry, so that = A 7 = 0. However, 
the scalar potential is now CP-violating unless In^Ag^m^] 2 ) = 0. In what follows, we shall assume 
that the latter condition is satisfied. In this case, one can rephase such that m ? 2 and A 5 are both 
real. That is, the CP-conserving scalar potential of interest is given by 

V = j$i + ml 2 $^$ 2 - m ? 2 [$|$ 2 + $J,$ 4 ] + AAi($J$i ) 2 + ±A 2 ($|$ 2 ) 2 + A 3 ($j$i)($3,$ 2 ) 


+A 4 ($i$ 2 )($ 2 < I > i) + 5A5 [($i$ 2 ) 2 + ($ 2 $i 


(41) 


where all scalar potential parameters are real. We denote this basis of scalar fields as the Z 2 -basis. 

If A 5 < | rri ? 2 1 / ("^l ^ 2 ), then the vevs are also real in light of Eq. (0), which implies that the vacuum 
is also CP-conserving. Otherwise, there is the potential for spontaneous CP-violation. To guarantee 
that the vacuum is CP-conserving, we check the conditions given by Eq. ([22]) . Setting Ag = A 7 = 0 in 
Eq. (1231) yields 

Im(ZgZ 7 ) = AA 5 (Ai - A 2 )s|g sin 2 £ . (42) 

The other two conditions given in Eq. (1221) simplify considerably when Ag = A 7 = 0, 

I m (Z§Zq) = — \\§slp sin 2 £-| (Ai — A 3 — A 4 ) 2 c| — (A 2 — A 3 — A 4 ) 2 s^ + ^As(Ai — A 2 )sip cos 2 £ — A 2 c 2 / g|> , 

(43) 

Im(Zg Z‘l) = — \\§slp sin 2 £<| (A 4 — A 3 — A 4 ) 2 s^ — (A 2 — A 3 — A 4 ) 2 c^ + ^As(A 4 — A 2 )s? ( g cos 2 £ + A?c 2( g , 

(44) 

where cos£ = m? 2 /(Au 4 7 ; 2 ) if A 5 > |m? 2 |/('Wiu 2 ) and sin£ = 0 otherwise. 

Although the vacuum is CP-violating for generic values of the softly broken Z 2 -invariant 2HDM 
scalar potential when sin£ 7 ^ 0, special cases can arise in which sin£ 7 ^ 0 and yet the vacuum is CP- 
conserving. We have already encountered one such example when m\ 2 = Ag = A 7 = 0 and A 5 > 0, in 
which case the scalar potential is minimized for £ = 77 t. More generally, consider the case of m\ 2 / 0, 
Ai = A 2 , A 5 = |A 4 — A 3 — A 4 | and Ag = A 7 = 0. In this case, Eqs. (1421) (1441) all vanish, which implies 
that a real Higgs basis exists. Indeed, for this particular example, Eqs. (fl 8 l) (1201) yield 

Z 2 

-y- =±\\ b slp(l ± cos 2£), Zg =-Z 7 = ±AA 5 s 2/3 [(1 ±cos2£)c 2/3 + isin2£] , (45) 
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for A5 = ±(Ai — A3 — A4). Thus, one can rephase the Higgs basis field H 2 to render Z 5 , Z 6 and Z 7 real. 
Henceforth, we assume that a Higgs basis has been chosen such that all the are real and sin/; = 0. 

The most general parametrization of the CP-conserving 2HDM scalar potential and vacuum is 
easily specified in the Higgs basis. To determine the constraints on the Higgs basis parameters that 
guarantee the existence of a softly-broken Z 2 -invariant, CP-conserving scalar potential, one simply 
sets Ag = A7 = sin £ = 0 in Eqs. (fTTTl (1201) . It is also convenient to redefine tan/3 as the ratio of vevs. 
That is, henceforth we shall define tan/3 = = e l ^v 2 /v\ = ±v 2 /v\ for £ = 0 and £ = 7 r, 

respectively. In particular, the Zi are now given by 

Z\ = AiCg + \ 2 Sfi + 5 A 345 s\p , 


Z 2 — AlS^g + A2 Cp + ^3455^ , 

z i = 14/3 [Ai + A 2 - 2A345] + A ? ;, (for i = 3,4 or 5 ) 
Zq = — \s2p\\i(?p — A2 s 2 p — A345 c 2 j a] , 

Z-j = — \s 2 p\\\s 2 p — \ 2 Cp + A345C2/3] , 


(46) 

(47) 

(48) 

(49) 

(50) 


where A 345 = A 3 + A 4 + A 5 . Since there are five nonzero Aj and seven nonzero Z % . there must be two 
relations. One can check that the following two identities are satisfied: 

Z 2 = Z\ + 2(.Zg + Z 7 ) cot 2/3 , (51) 

-^345 = Z\ + 2Zg cot 2/3 — (Zg — Zj) tan 2/3 , (52) 

where 

-^345 = Z 3 + Z4 + Z 5 . (53) 

Eliminating tan 2/3 yields the following relation among the Z t , 

(Zi — Z 2 ) [Z 1 Z 7 + Z 2 Zq — Z^(Zq + Z 7 )J + 2(Z 6 + Z 7 ) 2 (Zq — Z 7 ) = 0 . (54) 

Indeed, in a general CP-conserving 2HDM, if Eq. (154(1 is satisfied then there must exist a basis in 

which Ag = A 7 = 0. The corresponding value of tan/3 (which specifies the basis) can be determined 
from either Eq. m or (H 2 j). At this stage, the parameter tan/3 is a pseudo-invariant, since its sign 
can be flipped by redefining H 2 —>• — H 2 [e.g., 4>2 —>• — 4*2 and £2 —>• £2 + tt in Eq. ([H|>] . It is convenient 
to work in a convention where the ratio of vevs is non-negative, in which case we can take 

0</3<±vr. (55) 

In this convention, the signs of the pseudo-invariants Zq and Z 7 are now fixed, since we can no longer 
flip their signs by redefining H 2 —>• —H 2 . This implies that the relative sign of sp_ a and cp_ a is fixed 
by Eq. ([37]) . in which case the sign of cg_ a is determined by combining Eqs. (1371) and (l38l) . 

Finally, in the CP-conserving, softly-broken Z 2 -invariant 2HDM, it is convenient to introduce the 
squared-mass parameter, 

_ 2 _ 2 m^ _ m2 + ^2 _ (56) 

sm 2/3 

One can express m 2 in terms of Y 2 , Z\ and Zg, 

m 2 = Y 2 + \Z\v 2 + Zgu 2 cot 2/3 . (57) 

Then, combining Eq. (1571) with Eqs. (f26|) . (1551) and (1361) yields 

Z 7 v 2 = (m 2 h - m 2 H )sp- a cp- a + 2 cot 2/3 [ m 2 H s 2 p_ a + m 2 h c 2 p_ a - m 2 ] . (58) 
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2.4 Special forms for the Higgs—fermion Yukawa couplings 

We next turn to the Higgs-fermion Yukawa couplings. One starts out initially with a Lagrangian 
expressed in terms of the scalar doublet fields (i = 1,2) and the interaction-eigenstate quark fields. 
After electroweak symmetry breaking, one can identify the 3x3 quark mass matrices. By redefining 
the left and right-handed quark and lepton fields appropriately, the quark and charged lepton mass 
matrices are transformed into diagonal form, where the diagonal elements are real and non-negative. 
The resulting Higgs-fermion interaction Lagrangian in terms of the quark and lepton mass-eigenstate 
fields, U = ( u,c,t ), D = ( d,s,b ), N = (, zx T ), and E = (e,/r,r), is given by 

- Cy = U L ^a*h U a Ur - D L K^~h U a Ur + U L K^b% f Dr + Dl^X ] Dr 

+N L <S>+h^E R + E L <S>° a h^D R + h.c., (59) 

where K is the CKM quark mixing matrix, h u,D,L are 3x3 Yukawa coupling matrices and there is 
an implicit sum over a = 1,2. The diagonal quark and charged lepton mass matrices are given by 
Mr = {v\h[ + v^h^) /\/2, where F = U,D,E. However, the couplings of the neutral Higgs bosons to 
the fermions are not flavor-diagonal. Thus, Eq. (15911 would yield large tree-level Higgs-mediated flavor 
changing neutral currents (FCNCs) which is in conflict with observed data. 

In a general extended Higgs model, tree-level Higgs mediated FCNCs are absent if for some choice 
of basis of the scalar fields, at most one Higgs multiplet is responsible for providing mass for quarks or 
leptons of a given electric charge, as first pointed out by Glashow, Weinberg and Pascos (GWP) [351136] , 
This GWP condition can be imposed by a symmetry principle, which guarantees that the absence of 
tree-level Higgs-mediated FCNCs is natural. By an appropriate choice of symmetry transformation 
laws for the fermions and the Higgs scalars, the resulting Higgs-fermion Yukawa interactions take on 
the required form in a specific basis of scalar fields. The symmetry also restricts the form of the Higgs 
scalar potential in the same basis. These considerations were first applied in the 2HDM in Refs. m 
and [38] . 

The GWP condition can be implemented in four different ways [51H54] : 

1. Type-I Yukawa couplings: h\ r = hf = hf = 0, 

2. Type-II Yukawa couplings: h\ = hJJ = hlj = 0. 

3. Type-X Yukawa couplings: h\ = hf = = 0, 

4. Type-Y Yukawa couplings: h\ ] = h® = /if = 0. 

The four types of Yukawa couplings can be implemented by a discrete symmetry as shown in Table [T] 



$1 

4>2 

Ur 

Dr 

Dr 

Ul , E>l, Nr, E l 

Type I 

+ 

— 

- 

— 

- 

+ 

Type II (MSSM like) 

+ 

- 

— 

+ 

+ 

+ 

Type X (lepton specific) 

+ 

- 

— 

- 

+ 

+ 

Type Y (flipped) 

+ 

— 

— 

+ 

— 

+ 


Table 1: Four possible Z 2 charge assignments that forbid tree-level Higgs-mediated FCNC effects in 
the 2HDM [M]. 
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hUU 

hDD 

hEE 

HUU 

HDD 

HEE 

iAU-i^U 

iAD^/ 5 D 

iAE'y^E 

Type I 

cos a 
sin (3 

cos a 
sin (3 

cos a 
sin (3 

sin a 

sin (3 

sin a 
sin f3 

sin a 
sin (3 

— cot /3 

cot /3 

cot f3 

Type II 

cos a 
sin (3 

sincx 
cos (3 

sin a 
cos /3 

sin a 
sin /3 

cos a 
cos (3 

cos a 
cos (3 

— cot /3 

— tan/3 

— tan/3 

TypeX 

cos a 
sin (3 

cos a 
sin /3 

sin a 

cos /3 

sin a 
sin (3 

sin a 

sin (3 

cos a 
cos (3 

— cot /3 

cot /3 

— tan /3 

Type Y 

cos a 
sin (3 

sincx 
cos /3 

cos a 
sin (3 

sin a 
sin /3 

cos a 
cos (3 

sin a 
sin (3 

— cot /3 

— tan/3 

cot (3 


Table 2: Neutral Higgs-fermion couplings in the 2HDM subject to the Z 2 symmetries given in Table[TJ 
The couplings listed above are normalized relative to the SM Higgs couplings to UU, DD , and EE. 

The neutral Higgs Yukawa couplings (relative to the corresponding couplings of the SM Higgs boson) 
are conveniently summarized in Table [2] for the four possible implementations of the GWP condition. 

The imposition of the discrete symmetry also restricts the form of the Higgs scalar potential given 
in Eq. © by setting m\ 2 = \6 = A 7 = 0. As discussed in Section 12.31 the condition mf 2 = 0 can 
be relaxed. In the case of a softly-broken ^-invariant 2HDM with m\ 2 7 ^ 0, Higgs-mediated FCNCs 
are still absent at tree-level, although they are generated at one-loop order. Nevertheless, the size of 
these FCNCs may be phenomenologically acceptable, depending on the region of the 2HDM parameter 
space. This motivates us to focus on the Higgs scalar potential of the form given in Eq. (1411) . Note 
that the parameter tan/3 = is defined in terms of the vevs with respect to the Z 2 -basis of 

scalar fields, where the discrete t L 2j symmetry of the Higgs-fermion Yukawa interactions is manifest 
(i.e., where three of the six Higgs-fermion Yukawa matrices vanish). Indeed as previously advertised, 
the parameter tan j3 has been promoted to a physical parameter of the theory. 

2.5 Parametrizing the softly-broken Z 2 -symmetric CP-conserving 2HDM 

The scalar potential in the Z 2 -basis [Eq. (1411) ] is governed by three squared-mass parameters and five 
dimensionless quartic coupling parameters. Minimizing the scalar potential, we re-express rri 2 j and 
in terms of v = 246 GeV and tan/3. Excluding v (which determines the W and Z masses), we 
are left with seven real parameters: 777 2 2 , tan/3 and the A, (i = 1,2,... ,5). In the Higgs basis, the 
counting is also straightforward: after imposing the scalar potential minimum conditions the relevant 
real parameters (excluding v) are Y 2 and (i = 1,2,..., 7). Imposing one relation, Eq. (|54l) . to 
guarantee the existence of a basis where Xq = A 7 = 0 , we are again left with seven real parameters. 

A more physical choice of parameters would consist of a, /3 and the four scalar masses, nth, 
77177 , 777.^4 and m H ±. This leaves one additional parameter, which is usually chosen to be 777 2 2 or A 5 
(cf. Appendix D of [47]). In many of the previous studies of the 2HDM parameter space ]24fl26li29] - 
[!Hi33]j3i] , scans were performed over the parameters 777 / 1 , 77777 , 777 , 4 , "7_/7±, a and /3. Acceptable 
points in the scan must satisfy unitarity and perturbativity constraints [551156] , However, for random 
choices of m\ 2 and masses of the three non-SM-like Higgs bosons, one finds that the unitarity and 
perturbativity constraints are often violated. This is easily understood by examining the decoupling 
limit where Y 2 ^$> v. In this limit, 77777 ^ v in which case cp- a —> 0 [cf. Eq. (1391) ] and the properties of 
h approach that of the SM Higgs boson. Moreover, eqs. (ITU) . (1261) . (1361) and (1561) yield squared-mass 
differences, 


m\ ~ m 77 ± = 5(^4 — Z 5 )v 2 , m 2 H — m 2 A ~ Z 5 v 2 , 777 2 — 777)4 = ■ ( 60 ) 
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The conditions of unitarity and perturbativity constrain these squared-mass differences to be of 0(cu 2 ), 
where c 10 . 

To maximize the efficiency of costly scans over the 2HDM parameter space (e.g. in a global fit), 
one should define the input parameters of the model to include the mass of the observed SM-like Higgs 
boson and at most one additional Higgs mass. The other two Higgs masses and m 2 are determined 
in terms of the quartic Higgs self-coupling parameters. As long as the magnitudes of the quartic 
parameters are constrained to be less than 0(10), the scan over the other Higgs masses and m 2 will 
not violate the unitarity and perturbativity constraints. 

2.6 The hybrid basis of parameters 

In light of the discovery of a SM-like Higgs boson, we have information on two of the parameters of 
the 2HDM. Since the observed boson is not CP-odd, we shall identify it with either h or H. Given 
that the hVV and HVV couplings (V' = W ± or Z) relative to the SM Higgs boson Hsm are given by 

9hVV = 9h SM VV s l3-a , 9HVV = 9h SM VV c P~a > ( 61 ) 

the values of the SM-like Higgs boson mass (m^ or mu') and cp— a are already known with some 
precision. 

We therefore propose a “hybrid” strategy for specifying the input parameters for the softly-broken 
Z 2 -invariant 2HDM. The masses of the physical CP-even Higgs scalars are given directly as input 
parameters together with cg_ Q , which determines the phenomenologically important couplings of the 
CP-even scalars to the W± and Z bosons [cf. Eq. (1611) ]. and tan/3 which specifies the basis of scalar 
fields where the discrete symmetry of the Higgs-fermion Yukawa interactions is manifest. In addition, 
we specify the real Higgs basis self-coupling coefficients Z 4 , Z$, and Z-j as input parameters. We 
henceforth designate the input parameter set {mh, mu, cg_ Q , tan/3, Z 4 , Z 5 , Z 7 } as the hybrid basis of 
parameters as indicated in Tabic [3l In the convention adopted in Eq. (|38j) . we hav^l 

0 < S/ 3 _ a < 1 , and — 1 < cg_ Q < 1 , (62) 

in addition to the convention of Eq. (|55|) which implies that tan /3 is non-negative. 


Parameter 


m h 

m H 

C(3—a 

tan/3 
Z 4 , Z 5 ,Z 7 


Description 

Mass of the light CP-even Higgs boson 
Mass of the heavy CP-even Higgs boson 

cos(/3 — a), which determines the couplings of the CP-even Higgs bosons to VV 
Ratio of vacuum expectation values in the basis with manifest Z 2 symmetry 
Quartic couplings in the Higgs basis of 0(1) 


Table 3: 2HDM input parameters in the hybrid basis. By convention, we take tan /3 to be non-negative. 
For a more detailed description of their phenomenological relevance, see the text. 

4 Different convention choices for /3 — a appear in the 2HDM literature. For example, the convention 0 < cp- a < 1 
and — 1 < sp- a < 1 is employed internally by the 2HDMC code 14211431 . Another common choice in the literature is to 
take —^7r < a < ^ 7 r. Indeed, one is always free to define either a or /3 — a modulo 7r. It is a simple matter to translate 
among the various conventions. For any 2HDM parameter point (a,/3) where 0 < j3 < ^7r, compute the values of sp- a 
and cp- a . Then, to convert to the convention of Eq. (1(1121) . simply replace ( sp- a ,cp- a ] —¥ (-sp^ a ,—cp- a ) if sp- a is 
initially negative. Only the relative sign of sp- a and cp- a is physical under the independent convention of non-negative 
tan/I (the latter is universally employed in the 2HDM literature), as explained below Eq. (1551) . 
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With these seven input parameters, one may compute the real Higgs basis self-coupling coefficients 
Z\, Z 2 , Z 3 and Zq as well as the CP-odd Higgs and charged Higgs scalar masses and m H ± as 
follows. First, we employ Eqs. (fMl) and (1551) . 


_ s}_ a m 2 h + c 2 p_ a m 2 H 

O 1 

v z 

(m 2 - m 2 H )sfi_ a c^ a 

— - 9 - 

V z 

In order to ensure that unitarity and perturbativity constraints are not violated, one should choose 
the input parameters mn and cp- a such that m 2 H cp- a ~ O(ov 2 ), where c /$ 0 ( 10 ) in order to avoid 
a potentially large value of Zq. Next, we use Eq. ifTil) to compute Z 2 , 


(63) 

(64) 


C 2 


jn\ 


+ C 0-a m 


H 


+ 2 cot 2/3 



m _f/) S /3—oT/3—a _|_ £ 

V 2 ‘) 


(65) 


There is a danger that Z 2 may become too large if j3 is near 0 or - 77 r. This is not a true singularity since 
in the formal limit of tan 2 (3 = 0, it follows that Zq = Z 7 = 0 [cf. Eqs. ([46jl and (07])] corresponding 
to the inert 2HDM limit discussed briefly at the end of this Section. In practice, we restrict our scans 
over regions of tan f3 such that Z 2 is never too large. 

The parameter Z 3 is determined from Eq. (1521) . 


^3 = 


s}_ a m 2 h + c 2 p_ a m 2 H 


+ (2 cot 2/3 — tan 2/3) 


(' m\ - m 2 H )sp a c 0 c 


+ Zj tan 2/3 — Z± — Z 5 . ( 66 ) 


For this quantity, values of j3 near 0, \ir or appear to be problematical, yielding potentially large 
values of Z 3 . The case of (3 near 0 or ^ 7 r has already been mentioned below Eq. (1651) . The case of 
f3 = \ 7 T is not a singular limit, as it corresponds to Z\ = Z 2 and Z§ = Z 7 with Z 345 an independent 
parameter [cf. Eqs. (Hbl) ([50]) ]. In particular, for /3 = | 7 r, we must replace Z 7 with Zj, as an input 
parameter in the hybrid basis. We can sidestep this special parameter regime by avoiding values of 
tan (3 too close to 1. 

Observe that Eqs. (1M1) . (1651) and (1661) are consistent with the pseudo-invariant nature of tan/3, 
S/ 3 - a cp- a , Z§ and Z 7 . Since we have established a convention in which sp_ a and tan/3 are both non¬ 
negative [cf. Eqs. (1381) and (1551) ]. it follows that that full parameter space of the model requires the 
consideration of all possible non-negative values of tan /3 and and all possible sign combinations 
of C/ 3 _ a , Zq and Zj, subject to the constraint of Eq. (1371) . 

The masses of the CV -odd Higgs and charged Higgs scalar masses are determined from Eqs. (I24|) . 
(1261) and (]36l) 


m 2 A = rn 2 H s 2 p_ a + m 2 h c}_ a - Z 5 v 2 , (67) 

m 2 H ± =m\ - i(Z 4 - Z 5 )v 2 . ( 68 ) 

As a side note, one can also determine m 2 from Eq. (1581) . 

m 2 = m 2 H s 2 p_ a + rn 2 h c 2 p_ a + \ tan 2/3 (Z 6 - Z 7 )v 2 , (69) 

where Z e is given by Eq. (IM1) . 
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In practice, taking the parameters Z 4 and Z 5 as input parameters is not very intuitive. Clearly 
a more physical approach is to take tua and m H ± as input parameter and then compute Z 4 and Z 5 
using Eqs. (jS7D and (1U51) . Then these parameters can be employed in evaluating Z 3 via Eq. (1UU1) . 
The only danger with such an approach is that a poor choice of rriA and m^± will yield values of Z 4 
and Z§ that violate unitarity and perturbativity constraints, as discussed in Section 12.51 In practice, 
we first employ our hybrid basis to perform our parameter scans. Once we have identified suitable 
parameter regimes, we will then take a more physical approach by employing tua and mn± to define 
the corresponding benchmark scenarios. 

The parameter Z 7 is also not particularly intuitive, so one could advocate taking m 2 as the 
input parameter and then compute Z 7 using Eq. (1691) . However, m 2 is not a physical squared-mass 
parameter, so there is no particular advantage for adopting this strategy. A poor choice of m 2 would 
yield a value of Z 7 that violates unitarity and perturbativity constraints. Indeed, the CT’-conserving, 
softly-broken ^-invariant 2HDM always requires one extra parameter beyond the four physical Higgs 
masses, a and f3. There is an advantage to choosing this parameter to be dimensionless for the reasons 
noted above. Another possible choice for this parameter that is often found in the literature is A 5 . 
Using Eqs. (|48]) - (|50|) . one can check that 

A 5 = ^5 + 5(^6 — Z 7 ) tan 2/3 , (70) 

so one can always trade in A 5 for Z 7 and vice versa. 

The case of m 2 = 0 is special and corresponds to a 2HDM with an exact Z 2 symmetry. In this case, 
Zj is determined by Eq. (1691) in terms of m^, m# and cg_ Q . It is noteworthy that one cannot take mu 
arbitrarily large in the Z 2 -invariant 2HDM without violating unitarity and perturbativity bounds. 1 
Thus, there is no consistent decoupling limit where the heavy Higgs states H , A and H± decouple 
from the low-energy theory. In contrast, in the softly-broken ^-invariant 2HDM with m 2 7 ^ 0, the 
decoupling limit exists in which H , A and are all very heavy with squared-masses of 0(jn 2 ) 
[cf. Eq. (1571) ]. while keeping all Higgs self-coupling parameters bounded. 

Finally, the inert 2HDM |57W59j . corresponding to Z§ = Z 7 = 0 and Type-I Yukawa couplings EQ], 
must be treated separately. In the inert 2HDM, we have either (3 = 0 or (3 = and either cg_ a = 0 
(in which case h is identified as the SM-like Higgs boson) or sp_ a = 0 (in which case H is identified 
as the SM-like Higgs boson). In both cases, the SM-like Higgs boson is denoted by h-sM and the other 
inert neutral Higgs bosons will be denoted by Hi and Aj. The states Hi and A/ are relatively CV- odd, 
although (despite the notation) the individual CU-quantum numbers of these two states are not well 
defined. By convention, we shall define Hi to be the heavier of the two neutral inert Higgs scalars. 
Then, using Eqs. (|24fl and (1251) . the charged and neutral inert Higgs masses are given by 

m ‘ji± = Y2 + ^Z^v 2 , m ‘H I ,A I = ^2 + 5(^3 + ^4± \Z§\)v 2 . ( 71 ) 

The relevant input parameter set in this case is {m 2 gM = Ziv 2 , m 2 H ±, Z 2 , Z 3 , Z 4 , Z 5 }. Note that 
parameters Z 2 and Z 3 , which govern the self-couplings of the inert scalar doublet, must be provided 
separately since Eqs. fl65j) and (f 66 |) are not applicable in this case. We shall not consider the inert 
2HDM further here, but phenomenological constraints and benchmark scenarios for this model have 
been recently discussed in [60] . 

s In light of Eq. (I T 0 11 . the quantity (Zq — Zt)ta,n2l3 ~ 0(1) for all values of tan f) as long as the A i ~ 0(1). 
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cos(j0 — a) 


Figure 1: Example 2HDM parameter regions respecting perturbative unitarity and stability constraints 
(green) for mu = 300 GeV (left) and mu = 600 GeV (right), Z 4 = Z 5 = —2 and Z 7 = 0. 


3 Benchmark scenarios 

3.1 Numerical analysis 

To capture the 2HDM phenomenology that is interesting for upcoming LHC Higgs searches, we have 
performed numerical scans over the softly broken ^-symmetric 2HDM parameter space using the 
hybrid basis of parameters for inputs for the scans @ Based on these results, we identify a set of useful 
benchmark scenarios that we will now describe in more detail. Most of the scenarios are presented 
in the form of 2 -dimensional benchmark planes, whereas some of them are in the form of lines where 
only one parameter is varied. As a general rule, these scenarios could be made into benchmarks of 
higher dimensions by promoting additional parameters with fixed values to vary within certain ranges. 
Below we discuss a few examples of how this could be done. 

For the numerical evaluations we use the code 2HDMC [42fil3] (v. 1.7.0), where the hybrid basis of 
parameters has been implemented according to Section 12.61 Constraints on the quartic couplings from 
(absolute) vacuum stability and S-matrix unitarity are evaluated at the input scale. For the case with 
Z2 symmetry (as is implicit already in the definition of the hybrid basis), the condition that the Higgs 
potential is positive definite is equivalent to the well-known relations mi 

Ai > 0, A 2 > 0, A 3 > —A/A 1 A 2 , A 3 + A 4 — |A 5 1 > —y / AiA 2 . (72) 

For the unitarity constraints, we impose an upper limit which corresponds to saturation of the unitarity 
bound by the tree-level contribution. Using the formulation of [56], this is equivalent to the constraint 
|A| < 16 7 T on individual matrix elements. As an example, the regions in (cg_ Q , tan f3 ) compatible 
with theoretical constraints from perturbative unitarity and positivity of the potential are shown in 
Fig. [D for fixed values of the input parameters in the hybrid basis: mh = 125 GeV, mu = 300 GeV 
(left) and mu = 600 GeV (right), ^4 = ^5 = —2 and Z-j = 0. The theoretically allowed parameter 
space (shown in green) shrinks with increasing mu , owing to the existence of a proper decoupling limit 
where cp_ a is forced towards zero as mu 3> v [cf. eq. (1391) ]. As can be seen from this figure, there is 

6 The parameter scans presented in Section [3] have employed either Type-I or Type-II Yukawa couplings. Parameter 
scans with Type-X or Type-Y Yukawa couplings will be treated elsewhere. 
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an allowed (green) patch that survives for large values of cg_ Q even for mn = 600 GeV. This region 
is removed at higher masses by the unitarity constraint (which cuts in from the upper right corner). 
The constraints are invariant under the transformation (cg_ Q , tan/3) —>• (—c^_ a , l/tan/3), reflecting 
the fact that these points are related by an interchange of the Higgs fields in the Z 2 -basis, 4>i <—> <I? 2 , 
or equivalently taking H 2 —>• —H 2 and cp- a —>• —cp- a in a convention where sp- a is non-negative 
[cf. Eq. (1311) ]. We shall not consider values of tan f3 < 1 in the following. 

Branching ratios and LHC cross sections for Higgs production are evaluated using 2HDMC and 
SusHi [ 6 lj in accordance with the recommendations of [62]. We recommend that numerical values 
used in an experimental analysis are also calculated with HIGLU [63] and HDECAY [64], since such a 
comparison can provide a first estimate of the theoretical uncertainties in the treatment of (missing) 
higher-order corrections. These programs provide cross section predictions for the heavier 2HDM 
Higgs bosons, H and A, in the dominant gluon fusion and bb associated production modes. One of the 
benchmark scenarios discussed below (Scenario E) allows for the production of heavy charged Higgs 
bosons. The cross section for pp —>• tbH~/tbH + depends only on m^± and tan /?. A detailed numerical 
analysis of this process has recently been presented in [65], with numbers that are applicable to our 
Scenario E. 

Constraints from direct Higgs searches at LEP, the Tevatron, and the LHC are evaluated using 
HiggsBounds [ 66 U 68 ] (v. 4.2.0), which selects for each parameter point the most sensitive exclusion 
limit (at 95% C.L.). The compatibility of the 2HDM with the observed 125 GeV Higgs signal is 
calculated in terms of a y 2 value taking into account the full LHC run-I data using HiggsSignals [69] 
(v. 1.3.0). To determine the viable parameter regions in a benchmark scenario with two free parame¬ 
ters, we demand compatibility with the best-fit point, which usually is in very good agreement with 
the SM, within 2 <7 (Ay 2 = y 2 — y 2 lin < 6.18). The best-fit (minimal y 2 value) is reevaluated for each 
benchmark scenario. We make no quantitative statements about the relative degree of compatibility 
between different scenarios and the data, but it is checked explicitly that the best fit points obtained 
in all cases lie very close to the SM predictions. Fig. [2] demonstrates how these constraints work on 
the 2HDM parameter space under the assumption that the observed LHC signal corresponds to the 
lightest 2HDM Higgs boson, h. Demanding that the Higgs rates are within 2 a of the measurements 
restricts cg_ Q to be close to the SM limit (cp- a ~ 0) in the case of Type-II Yukawa couplings. For 
Type-I couplings, the deviation from the SM limit can be somewhat larger. The shape of the allowed 
region arises from the dependence of the production and decay rates on the mixing angles, most 
importantly the total width which in the Type-II model is dominated by the hbb coupling. 

One additional set of constraints to keep in mind when designing viable benchmark scenarios 
are the electroweak precision tests, where in particular the oblique T parameter m can receive 
sizable contributions due to a large mass splitting between the non-SM-like charged and neutral Higgs 
states [TIKE!. In particular, 


aT = 


647T 


|j r (m^±,m^o) + T(m 2 H± ,m 2 H o) - F(m\ 0 ,m 2 H0 ) 


+c}_ a [J 7 (m 2 H± ,m 2 h 0 ) - J 7 (m \ 0 , m 2 h0 ) + T ( m 2 A0 , m| 0 ) - T ( m 2 H± , m 2 H0 )] j + O (g r2 ), (73) 

where a = e 2 /( 47 r) with the electromagnetic coupling e and the weak couplings g and g' defined in 
the MS scheme evaluated at mz- The function T is defined by 


= \{m{ + m 2 2 ) - 


2 2 
m\m\ 

2 2 

mf — 7712 


7717 


In I —2 

mi 


(74) 
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Figure 2: Direct constraints from LHC Higgs searches on the parameter space for the 2HDM Type-I 
(top) and Type-II (bottom) with mjj = 300 GeV (left) and mjj = 600 GeV (right). In both cases 
nih = 125 GeV, Z 4 = Z 5 = — 2 and Z 7 = 0. The colors indicate compatibility with the observed 
Higgs signal at 1 a (green), 2 a (yellow) and 3a (blue). Exclusion bounds at 95% C.L. from the non¬ 
observation of the additional Higgs states are overlaid in gray. Note that the constraints from Fig. Q] 
have not been imposed here. 


with J-(m 2 ,m 2 ) = 0. Typically, the 0(g' 2 ) contribution to T is numerically small. Hence, to ensure 
that the 2HDM contribution to T is within the current bound, T < 0.2 at 95% C.L. m , one must 
be close to the custodial limit. In particular, in the decoupling limit of the 2HDM, Eq. (1731) reduces 
to [75] 


aT ~ 


(m 2 ± - m 2 A ){m 2 ± - m 2 H ) 


A8n 2 v 2 m 


l + O 


m 


H 


(75) 


from which one can see that the 2HDM contribution to the T parameter will be sufficiently small if 
i\ < 0(v 2 ) or m 2 H± —m 2 H < 0(v 2 ). Note that Eq. d75j) can be rewritten as [cf. Eq. ([60]) ]: 


either m 2 H ± — mZ 


aT ~ 


(Zj — Z^){Zi + Z^)v 2 
192ir 2 m 2 H 


l + O 


m 


H 


(76) 


Consequently, the choice Z\ = +Z 5 is special in the sense that the corresponding Higgs sector contri¬ 
butions to the T parameter vanish in the decoupling limit. Several of the benchmark scenarios that 
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we propose below satisfy this property. 

Finally, the general 2HDM (without supersymmetry) can also be constrained by various low-energy 
(flavor physics) processes (see, e.g., [7?]). Since our main objective is to define scenarios capturing 
interesting LHC phenomenology, we will not be concerned with the details of these constraints, nor will 
they be explicitly applied in our numerical results. In all the following, it should be noted that there 
exist generic lower bound on the charged Higgs boson mass in the Type-II 2HDM, m^± > 480 GeV 
at 95% C.L., from measurements of the BR(P —> X^y) [78| . 

3.2 Scenario A (non-alignment) 

Our first benchmark scenario has the “normal” interpretation of the 125 GeV signal as the lightest 
CP-even Higgs boson, h, with SM-like properties. The h is SM-like in the so-called alignment limit@ 
where cg_ Q — >• 0, in which case the hVV coupling approaches the corresponding SM value. On the 
other hand, to allow for some interesting phenomenology of the heavier CV even state, H, we define 
the scenario with a non-alignment (cg_ Q 7 ^ 0) as allowed by the present constraints [cf. Fig. [2]. The 
scenario focuses on searches for the heavier CP-even state, H, in SM final states (including the H —>• hh 
decay). The remaining two Higgs bosons, A and H± (which are kept mass-degenerate), are decoupled 
to a sufficient degree to create a small hierarchy 

rrih = 125 GeV < mn < mA = rn H ±. 

For mu > 150 GeV, this can be achieved by setting Z 4 = Z 5 = —2, which leads to values of m#± 
satisfying the b -A s'y constraint for Type-II models. The value of cp_ a is fixed close to the maximum 
allowed by the LHC Higgs constraints, cg_ a =0.1 for Type-I and cg_ Q = 0.01 for Type-II couplings. 
Consequently, we keep mu and tan /3 as free parameters. As shown in Fig. [3l these parameter choices 



200 300 400 500 600 700 800 900 1000 

mu (GeV) 



Figure 3: Parameter space of the non-aligned benchmark Scenario A with Type-I couplings, cg_ Q = 0.1 
(left) and Type-II Yukawa couplings, cg_ Q = 0.01 (right). The color coding is the same as in Fig. [2l 

7 In the alignment limit I20I24I34I79H81I . the SM-like Higgs boson is approximately aligned with the neutral component 
of the Higgs basis field H\. In particular, in the limit where cp- a -A 0, we have h ~ \/2R eH® — v , and ~ Z\v 2 . In 
light of Eq. (1391) . this limit can be realized if either Z§ —> 0 or if mn rrih- The latter is achieved in the decoupling 
limit, whereas the former can be achieved independently of the choice of the non-SM-like Higgs masses. 
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Figure 4: Cross sections for gluon fusion gg —>• H at y/s = 13 TeV in the 2HDM relative to the 
prediction for a SM Higgs boson with the same mass in the alignment limit cg_ a —>• 0 with Type-I 
(left) and Type-II Yukawa couplings (right). 


lead to a very good fit to the light Higgs signal rates over a large fraction of the (mu, tan/3) plane. 
An alternative to choosing a fixed value for cp— a is to realize this scenario with a value that decreases 
with mu to emulate decoupling [cf. Eq. (1391) ]. For example, in the case of Type-I couplings, we have 
performed scans with cp- a = 0.1 x (150 GeV/m#) 2 . 

We now look at the predictions for the production and decay modes of the heavy CP-even Higgs 
in this scenario. Fig. [4] shows contours of the cross section for the dominant gluon fusion production 
mode, gg H versus the free parameters (mu, tan/3). As can be seen from this figure, the cross 
section is maximized as low values of tan/3 due to the suppression of the Htt coupling by 1/tan /3. 
This holds for both Type-I and Type-II couplings due to the universality of the couplings to up-type 
quarks. In the Type-II case, we also see an enhancement of the cross section at high values of tan /3 
where the bottom loops become dominant, although in this scenario the parameter space regions with 
a significant enhancement is already disfavored by LHC constraints. In the Type-I case, the decoupling 
property of the H couplings to fermions sets in at high tan /3 , which can lead to arbitrarily small values 
of the cross section. This scenario can therefore provide a useful benchmark to gauge and interpret 
experimental progress in probing smaller and smaller values of a X BR in particular channels. In the 
Type-II model, at high tan/3, the cross section for associated production with b quarks, bb —> H , can 
become important since it scales as tan 2 /3, whereas for Type-I couplings this process is suppressed. 

Turning now to the decay modes, the case of Type-I Yukawa couplings is shown in Fig. [5] for 
tan/3 = 1.5 and 7. The Type-I couplings gives a very distinct pattern of dominant decay modes, 
where the vector bosons dominate below 250 GeV. Above this value, the cascade H —>• hh opens 
up and quickly becomes dominant, while there are still substantial contribution from the decays to 
WW and ZZ. At low tan/3 this changes above the tt threshold, but for higher values of tan/3 this 
mode (like all fermionic decays) is suppressed. Scenario A with Type-I Yukawa couplings is a suitable 
benchmark to replace the SM as the standard parametrization of 2£ and M searches at higher masses, 
where the rates can basically be adjusted to anything below the present limit. Note that arbitrarily 
high values of mu are not allowed in this scenario for all tan (3 values, owing to the restrictive unitarity 
and stability constraints when cg_ a is non-zero [cf. Fig. [T], 
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mu (GeV) mu (GeV) 

Figure 5: Branching ratios of the Heavy Higgs boson, H, in scenario A with Type-I couplings for 
c^_ Q = 0.1, tan/3 = 1.5 (left) and tan/3 = 7 (right). Colors: H —> W + W~ (blue, solid), H —> ZZ 
(red, solid), H —>• hh (green, solid), H —>• tt (gray, short dash), H —>• bb (black, long dash), H —>• rr 
(gray, long dash) and H ^ gg (black, short dash). 



mu (GeV) mu (GeV) 


Figure 6: Branching ratios of H in scenario A with Type-II couplings for cg_ a = 0.01, tan/3 = 1.5 
(left) and cg_ Q = 0.01, tan/3 = 7 (right). The color coding is the same as in Fig [5j 


With Type-II couplings, the restriction on mu from theoretical constraints is much less severe, 
since the LHC data already forces the scenario to be defined with cg_ Q much closer to the alignment 
limit. This also has the effect of suppressing bosonic decay modes, which can be seen in Fig. [6] where 
fermionic modes (bb, tt and tt) are clearly dominant over the full mass range. This scenario can be 
used for combinations of searchers utilizing different fermionic decay modes. It is also useful to gauge 
performance of e.g. rr searches in a setting where only a single CV-even Higgs boson is present at a 
given mass (in contrast to the MSSM case where typically H and A are nearly mass-degenerate). It 
also avoids other complications associated with SUSY interpretations, such as the non-holomorphic 
A b corrections to the bottom Yukawa coupling and/or the presence of SUSY decay modes [821183]. 


21 



















3.3 Scenario B (\ow-m H ) 

Scenario B corresponds to a “flipped” 2HDM benchmark scenario. In this scenario both h and H are 
light, but it is the heavier of the two which has mjj = 125 GeV and is SM-like. Since rrif, < mn, the 
lighter of the two CP-even Higgs bosons must have strongly suppressed couplings to vector bosons to 
be compatible with direct search limits which forces sp_ a —> 0 [cf. Eq. (I61l) ]l*l 

Similarly to Scenario A, the hybrid basis quartic parameters are chosen to decouple the other two 
Higgs states, A and H , to a sufficient degree not to affect the phenomenology (here: Z/± = Z§ = —5). 
This is to ensure compatibility with the strong indirect constraints (from flavor physics) on a light 
charged Higgs boson with mass of order mjj, in particular for Type-II Yukawa couplings. In the 
parameter space region 90 < < 120 GeV LHC constraints (from h —>• 66 , rr) apply, which leads 

to an upper limit on tan/3. This limit depends only weakly on cg_ a in the limit cp_ a <C 1. We focus 
here on the benchmark line with tan/3 = 1.5 as a function of rrih, although it is straightforward to 
generalize this analysis to a benchmark plane by varying tan (3. 



m h (GeV) m h (GeV) 

Figure 7: Allowed parameter regions for the lightest 2HDM Higgs boson in Scenario B with Type-I 
Yukawa couplings (left) and Type-II couplings (right). The colors indicate statistical compatibility 
with the 125 GeV signal at la (green), 2a (yellow) and 3a (blue). The gray (dashed) region is 
excluded at 95% C.L. by constraints from direct searches at LEP and the LHC. 

Since this scenario implements a completely different interpretation of the 125 GeV signal, the 
assumptions for the parameter space fit presented in Fig. [2] are not fulfilled. In Fig. [7] we therefore 
show the corresponding results with the assumption of mu — 125 GeV for a region of parameter space 
near the SM limit, Icg-^l —>• 1. From this figure it can also be seen that the region where H —>• hh is 
open is very hard to reconcile with LHC measurements, we will therefore not consider this possibility 
any further. 

As suitable benchmark scenarios we select the case of exact alignment, ca_ a = 1, with either Type- 
I or Type-II Yukawa couplings. As a third scenario, we choose a non-aligned value of cp- a = 0.9, 
which however is only relevant for Type-I couplings. The dominant decay modes of the light, non- 
SM-like, Higgs bosons are similar for all the three different scenarios, with BR(6, —>• 66) ~ 75-80% and 
BR(/i —>• t + t~) ~ 8% the two most interesting from the phenomenological point of view. 

8 In light of Eq. (l40l) . the alignment limit of sp- a —» 0 is achieved in the limit of Z§ —> 0, in which case the heavier of 
the two CP-even Higgs scalars is approximately aligned with the Higgs basis held, H ~ \/2Rei/i — v, and ran — Ziv 2 
[cf. footnote [7]- In contrast to Scenario A, the decoupling limit can never be realized in Scenario B since at least one of 
the Higgs states is lighter than the SM-like Higgs boson. 
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3.4 Scenario C (CP-overlap) 

In this work we have restricted ourselves to benchmarks for a 2HDM Higgs sector with CP-conservation, 
while a detailed analysis of the hybrid basis extended to the CP-violating case is postponed to future 
work. Nevertheless, we now consider a scenario where overlapping CP-odd and CP-even Higgs bosons 
simultaneously have mass close to 125 GeV [Mj. Since the CP-odd Higgs boson does not couple to 
vector bosons at tree level, there are surprisingly few channels where it is possible to distinguish this 
scenario from the case with a single light Higgs, h. The most important channel where the CP-odd 
contribution to the total rate could reach 0 ( 1 ) is through gluon (bb) fusion, followed by the decay 
h/A — >• t + t~ . We shall therefore analyze this process in more detail. 

Requiring that = nrih (= 125 GeV) in the hybrid basis of parameters, Eq. (1671) yields 

Z 5 = (77) 

Requiring in addition that m H ± = mu, Eq. (1551) implies that 


— —Z§ 


2 {m 2 H - m 2 )c 2 0 _ a 


(78) 


The remaining quartic parameter, Z 7 , does not enter in the mass determination. To maintain a SM- 
like h, we focus on the alignment limit and set cg_ Q = 0. In this case, Eqs. (1771) and (1751) reduce 
to 


C 5 


-z 4 


m 


H 



2 

h 


(79) 


We fix Z 7 , which has only minor impact on the phenomenology, such that m 2 2 = m 2 A s 0 c 0 , or equiv¬ 
alently A 5 = 0. Setting mu = 300 GeV, this leaves tan /3 as the only remaining free parameter in 
Scenario C. Other choices for mu (or even varying its value continuously) would lead to a benchmark 
plane generalization of Scenario C. However, varying mu has no impact on the properties of the 
overlapping CP states (h and A). 

Scenario C can be considered with both Type-I and Type-II Yukawa couplings. In Fig. [51 we show 
the production cross sections in the hadronic modes gg —>• h/A and bb —» h/A (and their sum) for 
Type-I couplings. In this case the CP-even production dominates, except at low values of tan/3 < 1.6. 
Since there is no tan /3 enhancement of the bottom Yukawa, the 66 -induced production remains small 
and can be neglected. Combining the (total) production cross sections (gg + bb) with the branching 
ratios BR (h/A —>• rr), we form the inclusive quantity 



a(pp —>■ h/A) x BR (h/A —»• rr) 
[a(pp ->• h/A) x BR(h/A ->• rr)] SM 


(80) 


for the rr rate relative to the SM. Predictions for R TT in Scenario C are shown in Fig. [9l The figure 
shows that the total rate (green curve) approaches the SM value (R T t = 1) in the limit of high tan (3 due 
to the decoupling property for the A contribution with Type-I couplings. In the opposite limit, the CP- 
odd contribution becomes increasingly important and dominates for tan/3 < 1.6. In effect, the total 
predicted rr rate significantly exceeds the SM prediction. Experimental constraints limit the maximal 
total rate, with the currently measured rates being p^/ LAS = 2.llg g [85] and = 0.34T1.09 [56] FI 

9 These results correspond to the “boosted” category of ATLAS and to the “0 jet” results of CMS, both of which are 
expected to be dominated by gluon fusion production. 
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Figure 8: Hadronic cross sections at yfs = 13 TeV for production of the light CP-even Higgs boson 
(long dashes), the CP-odd Higgs A (short dashes), and their sum (green, solid) in Scenario C with 
Type-I Yukawa couplings. 



Figure 9: The total rr rate (adding gg and bb production modes), relative to the SM, from h (long 
dashes), A (short dashes) and their sum (green, solid) in Scenario C with Type-I Yukawa couplings. 
Right: the respective fractions of the inclusive rr rate resulting from h (long dashes) and A (short 
dashes). 


To determine the fractions of the total rate which can be due to the CP-even and CP-odd 2HDM 
components, we show in the right panel of Fig. [9] their relative contributions R l TT /R^. Due to the 
decrease in the cross section, the CP-odd component can be seen to decrease monotonously with 
increasing tan/3. This scenario therefore provides a consistent model to parametrize an arbitrarily 
small CP-odd admixture in the 125 GeV signal. 

With Type-II Yukawa couplings, the production cross sections for the CP-odd state, A, have a 
different dependence on tan /3. The effect of this change is shown in Fig. [101 As can be seen from 
this figure the CP-odd contribution now becomes much more important at high tan (3 , both due to 
an increase in gluon fusion which has a minimum around tan /3 ~ 4 and due to a contribution from 
the bb —>• A process with a cross section increasing as tan 2 /?. In the case of Type-II couplings the 
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Figure 10: Hadronic cross sections at y/s = 13 TeV for production of the light CP-even Higgs boson 
(long dashes), the CP-odd Higgs A (short dashes), and their sum (green, solid) in Scenario C with 
Type-II Yukawa couplings. 



Figure 11: The total rr rate (adding gg and bb production modes), relative to the SM, from h (long 
dashes), A (short dashes) and their sum (green, solid) in Scenario C with Type-II Yukawa couplings. 
Right: the respective fractions of the inclusive rr rate resulting from h (long dashes) and A (short 
dashes). 


total cross section does not approach the SM value for any value of tan /3 . This also has consequences 
for the total rr rate, as can be seen in Fig. [TTJ In this scenario the minimum rate, obtained for 
tan/3 ~ 3.5, is R TT — 1.5. While 50% above the SM expectation, this is not experimentally excluded. 
Interestingly, the similar magnitude of the CP-even and CP-odd contributions leads to fractions that 
are nearly equal over a large range in tan/3 (Fig. [11] right). This scenario is therefore particularly 
suitable as a “best-case” test of the possibility to distinguish CP-properties of the 125 GeV Higgs 
boson in the rr channel. 
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3.5 Scenario D (Short cascade) 

This scenario is constructed with a SM-like h by fixing cp_ a to be zero (exact alignment). The mass 
hierarchy can be modified allow for either one—or both—of the decay modes H —> ZA or H —» iT =t i7 =F 
to be open, thus resulting in a “small cascade” of Higgs-to-Higgs decays0 These decays can typically 
be made dominant in the mass window 250 GeV < mu < 350 GeV (below tt threshold). A recent 
search for H — > ZA by the CMS Collaboration [23] already places interesting constraints on the 
parameter space of this type of scenario. Other modes that can be potentially of simultaneous interest 
is H —> hh and H —» A A (when A is very light). 

We present realizations of Scenario D for all the interesting cases below. We define these scenarios 
with a single free parameter by fixing tan j3 = 2, but this assumption could easily be relaxed. For 
simplicity, we choose two of the three non-SM-like Higgs masses to be equal. In the hybrid basis of 
parameters, these mass degeneracies can be implemented by an appropriate choice of Z 4 and Z5 as 
follows [cf. Eqs. (l67l) and ([68]) ]: 


and cg_ Q = 0 = 

1 

II 

(81) 

and C/ 3 - a = 0 = 

0 

II 

(82) 

mA = <= 

=> Z4 — z§. 

(83) 


In addition to the short cascade that we are investigating, additional “exotic” decay modes may 
be accessible, such as e.g. A —>■ (when is light). However, with the possibility of non¬ 

degeneracy between the heavy Higgs bosons (H and A in this case) there is in principle no guarantee 
for these other exotic decay modes to co-exist at any appreciable rate. We therefore choose to discuss 
only the H decays here, but there could be additional modes of similar type to exploit. 

Starting with the case of only a low rriA, this can be realized in the hybrid basis of parameters by 
choosing Z 4 = — 1, Z5 = 1. For mu close to 250 GeV, the decay H —>• A A can be open, with a rate 
that can be adjusted by varying Zy. Here we choose Z 7 = —1, which satisfies stability requirements 
in the whole mass ranged Fig. [T2l shows the branching ratios of the heavy CP-even Higgs boson, H, 
into the interesting final states. As this figures shows, there is nearly no difference in the decay pattern 
for the two Yukawa types (although, of course, the production cross sections are different). One thing 
that can be noteworthy is that the cascade decay, in this case H —>• ZA maintains an appreciable 
branching ratio also beyond the top threshold, with percent-level rates up to mu ~ 380 GeV. In 
Type-II models it is possible to suppress H tt further by going to higher tan f3 , but we find for the 
chosen set of parameters that this requires adjusting Zy to maintain positivity of the potential. 

An alternative mass hierarchy with mu± < mA = nr u can be arranged by setting Z 4 = 2 and 
Z5 = 0 while keeping the remaining parameters fixed. This generates the new possible decay modes 
H —>• W±HT (and H -A- H + H~ for very low m#±). Note that these low m#± might be in a region 
which is disallowed by flavor constraints. The decay branching ratios in this scenario are shown in 
Fig. El Similarly to the case with light A, the cascade decay H — >• is dominant over the 

range 250 < mu < 350 GeV, and with a sizable tail towards higher masses. In this case percent-level 
branching ratios are obtained for mu — 420 GeV. 

Finally, by setting Z\ = Z5 the hierarchy of masses becomes mA = rn>H ± < m H, with the light 
CP-odd and charged Higgs degenerate in mass. In this case all the different cascade modes for H can 
be open at the same time. The resulting branching ratios for the chosen scenario is shown in Fig. fl4l 

10 Such decay modes have also been considered recently in 187] . 

u In contrast, the opposite sign choice, Zj = 1 leads to problems with positivity of the scalar potential. 
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Figure 12: Branching ratios of H in Scenario D (with low tua) for tan j3 = 2 with Type-I (left) and 
Type-II Yukawa couplings. The colors show H —>• ZA (blue, solid), H —>• AA (blue, short dash), 
H —>• tt (gray, dash) and H —>• bb (black, long dash) and H — >• tt (gray, long dash). 



Figure 13: Branching ratios of H in Scenario D (with low m#±) for tan (3 = 2 with Type-I (left) and 
Type-II Yukawa couplings. The colors show H —> (red, solid), H —>• H + H~ (red, short dash), 

H —>• tt (gray, dash) and H —>• bb (black, long dash) and H —>• rr (gray, long dash). 


Note that the widths for the leading-order predictions for the decays involving one Higgs boson and 
one gauge bosons, e.g. H —> ZA, are proportional to gauge couplings (and therefore relatively fixed), 
whereas the modes involving triple-Higgs couplings, such as H —>• AA, are proportional to scalar 
couplings appearing in the Higgs potential. For example 1^1 

3haa = ~ v [(^3 + Z^ — Z§)cp- a — ZiSp-a\ . (84) 

The strength of such modes is therefore more sensitive to the specific choice of the 2HDM parameters. 
12 A list of the triple Higgs couplings expressed in terms of the Higgs basis quartic parameters can be found in [34] . 
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Figure 14: Branching ratios of H in Scenario D with both tua and mu± low (Z 4 = Z 5 ) for tan/3 = 2 
with Type-I (left) and Type-II (right) Yukawa couplings. The color coding is the same as used in 
Figs. [T2l and fl3l 


3.6 Scenario E (Long cascade) 

The short cascade (Scenario D) was defined with a one-step decay involving two (or more) Higgs 
bosons. In Scenario E, we extend this to a “long cascade” where two-step decays involving all three 
non-SM-like Higgs bosons are possible, assuming that H is always the lighter of the three. Assuming 
first (El) that H± is heavier than both A and H (realised in the hybrid basis by choosing Z 4 = — 6 , 
Z 5 = —2 (with Zj always zero in this scenario), can give rise to a “long” cascade and a second 
complementary direct decay that is always present [55] . 


H* -tW^A-tW^ZH 

-> W ± H. 


(85) 


The other hierarchy, with 1714 > m H ± > mu is achieved by setting Z 4 = 1, Z 5 
to long cascades which are inverted compared to Eq. m m- 

A ^ H T ^W ± W T H 

A -> ZH. 


—3 (E2). This leads 


( 86 ) 


Interesting signatures in this scenario include a heavy Higgs boson H in the t + t~ or bb channels 
augmented with multiple leptons and/or additional jets. 

As in scenario D, the recent search for A —>• ZH [23] places constraints on the parameter space. It 
should be interesting to analyze the interplay between this mode and other available signatures in this 
scenario in more detail. In Tabled] we present the branching ratios for the various cascade decay modes 
in the two different incarnations of Scenario E for input masses mu = 200 GeV and mu = 300 GeV 
with Type-I Yukawa couplings. Note that for this low value of tan /3 there is no significant change 
in any of these branching ratios in a Type-II setting. The main difference between the two types are 
instead in the decay of H at the end of the cascade, where the final states with down-type fermions 
( b, t) are more favored with Type-II couplings. 


28 








Scenario 

Masses (GeV) 
uih tua rriH± 

Branching ratios 

H ± ^W ± A H ± ^W ± H A^ZH A^W ± H^ 

El.l 

200 402 532 
300 460 577 

0.053 0.79 0.62 

0.041 0.74 0.39 

El.2 

200 471 317 

300 521 388 

0.27 0.56 0.25 

0.026 0.50 0.20 


Table 4: Mass spectrum and Type-I Higgs branching ratios to bosonic decay modes in Scenario E. 


There are several things of generic interest that can be noted from Tabled) The heaviest (parent) 
Higgs boson has a rather high probability, 70-85%, to decay into a lighter Higgs and a vector boson. 
However, most of this decay goes directly into the lighter of the two states (H) due to the larger 
available phase space. The aggregated branching ratio for a “long” cascade is thus suppressed, reaching 
typical values at the (few) percent-level (up to 5-6% in some cases). Nevertheless, to be able to study 
these types of final states would be an intriguing possibility, that would contain a lot of information 
about the structure of an extended Higgs sector. 

3.7 Scenario F (Flipped Yukawa) 

The flipped Yukawa scenario is characterized by SM-like couplings for the light Higgs, h, except for 
the couplings to down-type fermions which has a change of sign relative to the SM. This scenario is 
realized with Type-II Yukawa couplings for values of (cg_ Q , tan/3) solving the equation [89] 

= —-—- = sin(/3 — a) — tan /3 cos(/3 — a) = — 1. (87) 

Shdd cos P 

This unconventional solution is what generates the second branch of 2 cr-allowed parameter space in 
Fig- E) (right), with valid solution for not too small values of tan/3. We define Scenario F by simply 
fixing the value of cq_ q from the chosen tan /3 using this relation, while the remaining parameters are 
given the same values as in Scenario A: = Z§ = — 2, Z-? = 0. 

The (small) effects on the h phenomenology induced by the opposite sign bottom Yukawa coupling 
have been studied in detail in [89]. However, the comparably large allowed values of cp- a also open 
possibilities for the decay modes of the heavier Higgs boson beyond what can be achieved near the 
alignment limit as considered, for example, in Scenario A. To illustrate this point, we show in Fig. [15] 
the decay branching ratios of H as function of mn and tan/3. As this figure shows, the decays into 
vector bosons can be sizable over a large part of the remaining parameter space, clearly beyond what 
is possible in the alignment case. This observation holds even in the region above the top threshold, 
which makes the Flipped Yukawa scenario a suitable benchmark to replace the SM in interpretations 
of future heavy Higgs searches in the dilepton and M final states over a larger mass range. 

3.8 Scenario G (MSSM-like) 

This scenario is inspired by the Higgs potential of the MSSM. The tree-level MSSM is defined by the 
following values for the quartic couplings (in the “SUSY-basis” where supersymmetry is manifest), 

Ai = A 2 = \{g 2 + g' 2 ), A 3 = \{g 2 - g 2 ), A 4 = -\g 2 , A 5 = A 6 = A 7 = 0 . ( 88 ) 

Following Eq. (1561) . the remaining parameter is given by m\ 2 = sin 2/3 in terms of the more usual 
MSSM inputs: the CV- odd Higgs mass, and tan/3. In the tree-level MSSM, the predicted mass 
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Figure 15: Branching ratios of H in Scenario F. The panels show H —» WW (upper left), H —>• ZZ 
(upper right), H —>• tt (lower left) and H —>■ hh (lower right). Contours for each branching ratio are 
indicated with colors (see legend). Gray regions are excluded at 95% C.L. from direct searches. 


of the lightest Higgs boson, rrih < mz, is not compatible with the LHC measurements. However, the 
CP-even mass matrix receives corrections beyond leading order. To a good approximation, the leading 
radiative corrections can be parametrized as an additional contribution to A 2 |90j , corresponding to 
the shift A 2 —>• A 2 + 5, while (sub-leading) MSSM contributions to the other quartic couplings are 
neglected. This is reminiscent of the approach pursued in m- Since the leading radiative corrections 
in the MSSM can most easily be specified in terms of the SUSY-basis, we use this instead of hybrid 
basis input for Scenario G. Working in this approximation, we define our MSSM-inspired 2HDM 
scenario by three parameters: m^, tua and tan /?. Using this input, together with the MSSM relations 
for the A, (i % 2) given by Eq. (jSHj) . we solve for the value of 5 necessary to reproduce the desired value 
of rrih for the chosen values of tua and tan (3. In practice, this is done with 2HDMC using an iterative 
procedure^! To satisfy the constraints of unitarity and perturbativity, we impose the additional 
condition A 2 < 47 t. 

Although this scenario is inspired by the MSSM, it is (like all our other benchmarks) completely 
defined within the CP-conserving, softly-broken ^-symmetric 2HDM. There is therefore no principal 
restriction to the Type-II structure for the Yukawa couplings, although we choose to stick to this 
familiar pattern here to keep the connection to the MSSM tree-level structure. 

13 This procedure can be accessed through the set_hMSSM method or by running the CalcHMSSM example program. 
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Figure 16: Allowed parameter space by direct Higgs search constraints in the “MSSM-like” 2HDM 
with Type-II Yukawa couplings. The color coding is the same as in Fig. [2j 


Using a fixed value mh = 125 GeV, we scan over the remaining parameters rriA and tan /3 to 
determine the viable regions of parameter space. The results are shown in Fig. [I6l where green color 
indicates regions compatible with mh = 125 GeV (and A 2 < 47t). The shaded (gray) regions show 
the excluded regions at 95% C.L. from direct Higgs searches. In particular the limit H/A —>• rr 
plays a very important role to constrain this scenario both for high and low values of tan /3 for values 
mA <2 mt- Constraints from heavy Higgs searches benefit in this scenario from the near-degeneracy 
of the heavy Higgs bosons H and A. 

The resulting phenomenology of this scenario is very similar to the MSSM Higgs sector in the 
absence of additional low-energy degrees of freedom, which in the MSSM could provide additional 
decay channels for the heavy non-SM-like Higgs bosons. This can be seen, for example, from the 
region allowed by the LHC Higgs measurements in Fig. [T6j In the absence of large radiative corrections 
affecting the CP-even Higgs mixing EB, the LHC Higgs data forces the MSSM into the decoupling 
limit. Hence, we find mA > 360 GeV at 95% C.L. almost independently of tan/3. A related analysis 
in the MSSM context at low values of tan/3 with a heavy supersymmetric spectrum that makes use 
of the framework described above has recently been presented in [92j . 

3.9 Additional scenarios for consideration 

There are a number of additional scenarios worthy of consideration that we have not included in this 
work. For example, one can consider the inert 2HDM with Type-I Yukawa couplings, which provides 
a plausible dark matter candidate [59 (.l93Lf96] . As noted at the end of Section 12.61 phenomenological 
constraints and benchmark scenarios for this model have been recently discussed in [60]. 

It is also possible to consider simplified models of 2HDM phenomena, where h is SM-like and one of 
the non-SM-like Higgs states among H, A and is significantly lighter than the two heaviest states 
(keeping in mind that the split spectrum must be consistent with constraints due to the T parameter). 
In this case, one can focus on the phenomenology of a single non-SM-like scalar. Scenarios A and B 
already provide examples of this type, in which A and H ^ are significantly heavier than h and H. But 
other scenarios could be considered that would feature the production and decay of A or alone. In 
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these latter scenarios, it would be appropriate to slightly modify the hybrid basis, which specifies two 
scalar masses as input parameters. When focusing on either A or H±, it would be more appropriate 
to specify the SM-like Higgs mass and one other scalar mass [m^ or m#±] as the two input scalar 
masses. Of course, one could then use Eqs. (1671) and (| 68 l) to determine the masses of the remaining 
non-SM-like Higgs scalars. 

4 Conclusions 

We have introduced a new “hybrid” basis to define input parameters in the general, CR-conserving, 
two-Higgs-doublet Model with a softly broken Z 2 symmetry in a way that naturally accommodates 
constraints from perturbativity and unitarity. In this basis, the input parameters are the two CP-even 
Higgs masses, rrih and m#, three quartic couplings defined in the Higgs basis, Z 4 , Z 5 , and Z 7 , and 
the mixing angles cos (/3 — a ) and tan /3 which are very relevant for the phenomenology. In this hybrid 
approach, tan /3 can be interpreted as specifying the basis where the soft ./^-breaking is manifest. 

The hybrid basis of parameters has been implemented as an input option for the public computer 
code 2HDMC. Using this setup, we have performed a scan and a detailed numerical analysis of the 
parameter space constraints and predictions of cross sections and branching ratios for different 2HDM 
configurations. Using the results of this analysis, we have defined a set of 2HDM benchmark scenarios 
that we deem relevant for the design and interpretation of Higgs searches beyond the SM at LHC 
run-II and beyond. The scenarios are defined to cover different aspects of the 2HDM phenomenology, 
including interpretation of the 125 GeV Higgs signal as either the light or the heavy CV-even Higgs 
boson. We also devise a new scenario with overlapping CP-even and CP-odd Higgses at 125 GeV, which 
could potentially be very useful for the interpretation of CP-studies in the rr final state, as well as a 
scenario where the Yukawa coupling of the Higgs to down-type fermions has the “flipped” sign (relative 
to the coupling to vector bosons). Another two sets of scenarios provide Higgs cascade decays as a 
prominent feature, giving rise to completely non-standard final states that may contain additional 
leptons and/or jets. Finally, we define the 2HDM equivalent of the decoupled MSSM, where the 
Higgs spectrum follows MSSM relations with the dominant radiative corrections to the Higgs masses 
completely captured in a shift of the quartic coupling A 2 appearing in the Higgs potential. Results 
presented for this scenario should be immediately familiar to anyone following MSSM Higgs searches 
at the LHC, which could make it useful for comparison and communication of results. 

We hope that the presented scenarios will be useful and serve to inspire new discoveries! 
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A. Parameter values for 2HDM benchmark lines and planes 


For convenience, this appendix contains a list of the input parameters that can be used to realize the 
2HDM benchmark scenarios described in this paper. 



m h (GeV) 

Scenario A (Non-alignment) 

m H (GeV) cg Q Z 4 

^5 

Z 7 

tan /3 

Type 

Al.l 

125 

150...600 

0.1 

-2 

-2 

0 

1 . 

..50 

I 

A1.2 

125 

150...600 

0.1 x ( 150GeV f 

V m H ) 

-2 

-2 

0 

1 . 

. .50 

I 

A2.1 

125 

150.. . 600 

0.01 

-2 

-2 

0 

1 . 

..50 

II 

A2.2 

125 

150.. . 600 

0.01 X t 1 50GeV\ 2 

V m H J 

-2 

-2 

0 

1 . 

..50 

II 


Scenario B (Low-m//) 



m h 

(GeV) 

rriH (GeV) 

C/3 —a 

^4 

^5 

Z 7 

tan (3 

Type 

Bl.l 

65. 

.. 120 

125 

1.0 

-5 

-5 

0 

1.5 

I 

B1.2 

80. 

.. 120 

125 

0.9 

-5 

-5 

0 

1.5 

I 

B2 

65. 

.. 120 

125 

1.0 

-5 

-5 

0 

1.5 

II 



m h 

Scenario C (CP-overlap) 

mu mA m H ± c p~a A 5 tan (3 

Type 

Cl 

125 

300 125 300 

0 0 1...10 

I 

C 2 

125 

300 125 300 

0 0 1...10 

II 



m h (GeV) 

Scenario D (Short cascade) 

m H (GeV) cg Q Z 4 Z 5 

Z 7 

tan /3 

Type 

Dl.l 

125 

250...500 

0 

-1 

1 

-1 

2 

I 

D1.2 

125 

250...500 

0 

2 

0 

-1 

2 

I 

D1.3 

125 

250...500 

0 

1 

1 

-1 

2 

I 

D2.1 

125 

250...500 

0 

-1 

1 

-1 

2 

II 

D2.2 

125 

250...500 

0 

2 

0 

-1 

2 

II 

D2.3 

125 

250...500 

0 

1 

1 

-1 

2 

II 


Scenario E (Long cascade) 



m h (GeV) 

m H 

(GeV) 

C/3 —Q 

z 4 

^5 

Z 7 

tan /3 

Type 

El.l 

125 

200 

...300 

0 

-6 

-2 

0 

2 

I 

El.2 

125 

200 

...300 

0 

1 

-3 

0 

2 

I 

E2.1 

125 

200 

...300 

0 

-6 

-2 

0 

2 

II 

E2.2 

125 

200 

...300 

0 

1 

-3 

0 

2 

II 


m h (GeV) 

Scenario F (Flipped Yukawa) 

m H (GeV) cg Q Z 4 Z 5 Z 7 

tan /3 

Type 

F2 125 

150... 600 sin 2/3 -2 -2 0 

5... 50 

II 


Scenario G (MSSM-like) 

rrih (GeV) m,A (GeV) tan/3 Type 

G2 125 90...1000 1... 60 H 

A i (i / 2) given by their MSSM values [cf. Eq. ([ 88 ])]: A 2 determined by rrih 
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